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INTRODUCTION 

Because of the extensive use of Electro-explosive de- 
vices (EED's) in Apollo subsystems, there is a need for 
methods of man rating these devices. One such method is the 
"Bruceton Method," developed at Bruceton, Pennsylvania, in 
1944, by the Statistical Research Group (1), Princeton Uni- 
versity. Since that time, there have been many extensions 
and modifications of the original Bruceton Method. While a 
few of these changes were for the better, the overall effect 
was the production of a number of contradictory documents. 

As a consequence, many experimenters have misused_or misap- 
plied the Bruceton Method; therefore, the results they ob- 
tained were invalid. 

The purpose of this report is to provide a reliable up- 
to-date guide for the application of the Bruceton Method to 
EED’s. A number of areas that have caused confusion in the 
past are discussed, for example: the concept of the all- 

fire and no-fire point, the response curve and what it means, 
the sample size required, the concept of confidence limits,....’ 
and the critical assumption of normality. 
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SYMBOLS 


u 

a 

m 

s 

A.P. 


A.P.C. 


N.P. 

N.P.C. 

S m 


the mean of the normal distribution 

the standard deviation of the normal distribution 

an estimate of u 

an estimate of a 

estimate of the All— Fire point 

estimate of the All-Fire point with 90$ confidence 
estimate of the No— Fire point 

estimate of the No-Fire point with 90$ confidence 
estimate of standard deviation of y 
estimate of the standard deviation of a 
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RESPONSE CURVE 


The purpose of the Bruceton test is to estimate a re- 
sponse curve. The response curve represents the percent of 
units that should fire at a given current. Typical response 
curves are shown in GRAPHS 1-4. The Bruceton test requires 
that the response curve be a cumulative normal curve. This 
curve, shown in GRAPH 5, is characterized by two parameter 
y , the mean, and cr , the standard deviation; i.e., if y 
and o are known the entire curve can be generated. Hence, 
the purpose of the Bruceton test can be stated more explicitly: 
the purpose of the Bruceton test is to obtain an estimate of 
y , namely jn ; and an estimate of a , namely s ; and 
with these estimates generate an estimated response curve. 

All-Fire and No-Fire Points 


The term all-fire point has come to mean that current 
at which 99.9 percent of the units will fire and the no-fire 
point has come to mean that point at which .1 percent of the 
will fire. Once estimates of y and o are obtained, 
an estimated response curve can be generated; thus, these 
two points can be located. For example, in GRAPH 5 it is 
seen that the all-fire point is 1.04 and the no-fire point 
is 3.56. The above definition of all-fire and no-fire points 
should be used regardless whether the experimenter is inter- 
ested only in the no-fire point and not the all-fire point 
and vice-versa. In this connection, this paper will not en- 
gage in semantics in trying to define a "success 1 or a "fail- 
ure"; th.e.s.e two words have been left out of the text due to 
the confusion they have caused in the past. If the experi- 
menter is consistent in his definition of all-fire and no- 
fire, there can oe no confusion. 
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GRAPH 3 

TtPJCAU RESPONSE CURVE FOR E L E C T R O - E X P U O S I V E DEVICES 
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Confidence Limits 

Since m and s are only estimates of u ana a 
respectively, the estimated response curve will differ from 
the true, but unknown, response curve. Hence, the estimates 
of the AP point and the NP point may be in error. In order 
to overcome these errors, the idea of confidence limits has 
been employed. This is explained best by an example. Suppose 
that in a metal container there were 10,000 white and black 
balls and that someone was interested in estimating the per- 
centage of white balls. He could count the number of whites 
in the container, but that would be rather tedious. Rather 
he takes a sample of 300 and counts the number of white balls 
in his sample; this number divided by 300 represents an esti- 
mate of the true proportion p of white balls in the con- 
tainer. Note that it is only an estimate of the true value. 
It— may be a very poor estimate if the sampling plan had de- 
fects in it. Now suppose he observed 235 white balls in his 
sample of 300. Then p , the estimate of p would be equal 
to '300' = But t * le ex Perimenter is not content with 

knowing that p = .783, he wants some confidence in this 
number. Depending upon what he wants, he can choose between 
three types of confidence limits. 

They are : 

(1) two-sided; i.e., he wants to know that with 100a 
percent confidence that the true proportion p 
lies between two numbers, say a; and a 2 


9 


4 


(2) one-sided to the left; i.e., he wants to know that 
with 100a percent confidence that the true propor- 
tion p is less than some number b . 

(3.) — one-sided to the right; i.e. , he wants to know that 
with 100a percent confidence that the true propor- 
tion p is greater than some number ci . 

Note that a is some number between 0 and 1 that represents 
how confident the experimenter wishes to be. Suppose that he 
choose a = .90. The three confidence limits would then be 
as follows: _ 

The experimenter would be 90 percent confident that the true 
value of the proportion p would 


(1) 

lie between .740 

and .822 i.e. , 

(2) 

b less than: i. 

e. , p <_ .817 

(3) 

be greater than: 

i.e. , p <_ 749 


The interpretation given to the above confidence limits is 
as follows: If the experimenter repeated his experiment 100 

times and constructed his confidence limits, the true proper 
tion p would lie within the prescribed limits 90 percent 
of the time. He would erroneously assume p is within the 
prescribed limits about 10 percent of the time. 

When confidence limits are used for the all- fire and 
points, tt is obvious that the most useful ones are 
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the one-sided confidence limits. The experimenter does not 
want to know that the all-fire point lies between ai and 
a 2 . He usually has some maximum current that can be supplied 
and he wants to make sure the all-fire point is well below 
this level. The same logic holds true for the no-fire point. 
Many of the documents followed in the past explain only the 
two-sided confidence limits. This results in their use re- 
gardless of. the situation. 

SAMPLE SIZE 

The number of specimens in the sample determines to a 
large extent the amount of information which can be derived. 

A limited amount of information can be obtained from small 
samples ; precise conclusions require large samples. If one 
wishes to locate a mean only within wide limits a small 
sample will suffice. If, however, it is .necessary to esti- 
mate the mean accurately, a large sample will be required. 

The fact that sensitivity data do not consist of exact 
measurements reduces the amount of information contained in 
them. Roughly speaking, the amount of information is at 
best about half that which would be available in measured 
data. That is, a sample of 100 explosions and non-explosions 
will contain information equivalent to that in a sample of, 
at most, 50 measured critical currents. Sample sizes will 
therefore have to be larger than in statistical investigations 
where exact measurements may be made. 

The method of statistical analysis of the data has been 
simplified considerably by the use of approximations which 
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become., poor for small samples. This puts a lower limit on 
the sample size if the present method of analysis is to be 
used. 


In the original document on the Bruceton Method (1), 
the authors state repeatedly that their method is not valid 
for small samples, especially when the sample size is less 
than 100. 


ixiS next section which contains the Bruceton procedure 
has been taken for the most part from the original Bruceton 
report (1). However, a number of changes have been made in 
order to up-date the original report. 

BRUCETON PROCEDURE 


Preliminary Remark 

A considerable amount of statistical investigation deals 
with quantities which are continuous variables, but which 
cannot be measured in practice. For example, in testing the 
sensitivity of electro— explosive devices to current, a com- 
mon procedure is to test a number of EED*s at a given current 
There are currents at which some devices will fire and others 
will not, and it is assumed that those which did not fire 
would have fired had the current been sufficiently large. 

The assumption is, therefore, that different devices have 
different "critical" currents, and that if a device is tested 
at a given current it will or will not fire depending on 
whether the current at which it was fired is greater than or 
less than the critical current of that specimen. Further- 
more, it is not possible to "creep up" on the critical current 
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of a specimen by testing it at successively greater currents 
until it fires because this procedure alters the character- 
istics of the explosive mixture and thereby changes its crit- 
ical current. 

A current is a simple measurable variable but it is 
clear that the critical current of electro-explosive devices 
are not directly measurable, in practice. All that one can 
do with a single device is determine whether its critical 
current is greater than or less than a chosen current at 
which the device is tested. 


Methods will be provided for estimating parameters (mean, 
standard deviation, percentage points) and for testing statis- 
tical hypotheses. The methods will not be valid for small 
samples as numerous approximations will be made which rapidly 
become less accurate as the sample size decreases from 100. 

The methods are based on the assumption that the critical 
currents (or some known function of them) have, a normal prob- 
ability distribution. It is important that this assumption 
be' reasonably well satisfied. 


If the currents are not normally distributed, it is 
desirable to transform them so that they will have a normal 
distribution. In biological research, it is usually found 
that the logarithms of the critical concentrations themselves 
are normally distributed. There is some evidence in explo- 
sives research that logarithms of critical currents are more 
nearly normally distributed than the critical cuj rents them- 
selves, but the evidence is not so conclusive or so gene.-.’ , : v 
accepted in this field. This is a question, therefore, which 
will usually require some preliminary investigation. Any 
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effort devoted to finding whether the critical currents are 
normally distributed, and if not, what function of the crit- 
ical currents is normally distributed., will be well worth- 
while. 


The Experimental Procedure 

The procedure, described in terms of the explosives 
experiment, consists of the following steps: 

(1) Choose a current h at which the first specimen 
will be tested, and an interval d which will be 
the distance between testing currents. 

(2) If the first specimen explodes when tested at cur- 
rent h , the second specimen will be tested at 
current h-d . If the first specimen does not ex- 
plode, the second specimen will be tested at cur- 
rent h+d . 

( 3 ) In general a specimen will be tested at a current 

d below the current at which the previous specimen 
was tested if that specimen exploded, and d aoove 
the current at which the previous one was tested 
if it did not explode. 

In this manner one will obtain a sequence of fires and 

no-fires, where the x's denote fires and 0's denote no- 
fires. 


14 



t. 


t 


Height 



Here the first specimen did no.t_.explode , so the second one 
was tested at h+d ; the second did explode, so the third 
was tested at h , the level just below h+d ; the third 
exploded, so the fourth was tested at h-d ; the fourth 
exploded, so the fifth was tested at h-2d ; the fifth did 
not explode, so the sixth was tested at h-d , the level 
just above h-2d. 


Choice of h and d 

On the basis of past experience with other explosives, 
it is usually possible to make rough estimates of the mean 
and standard deviation for the explosive to be tested. 3uch 
rough estimates will be useful in choosing h and d for 
the experiment to be performed. 
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If m and s are the estimated mean and standard devi- 
ation respectively, then good choices for h and d are 
simply h = m and d = s . These are not the best choices 
for all circumstances, but they will serve very well for most 
purposes. Any great improvement in them would require that 
the preliminary estimates be very accurate. 

Computation Estimates of the Mean and Variance . 

It is assumed that there is known or has been found by 
preliminary investigation a variable (the critical current, 
or some function of the critical current), which is normally 
distributed. Any reference to critical current in the fol- 
lowing discussion will be to these "normalized" critical cur- 
rents. The original testing current h and the interval d 
will be measured in terms of the normalized current. If, for 
example, the normalized current is the logarithm of the actual 
currents and the normalized interval is one, then corresponding 
to normalized current 0, 1, 2, 3, would be actual testing 
currents of 0, 1, 10, 100, 1000 amperes. 

Figure 6 is a chart showing the result of an actual test 
of 100 specimens of an explosive; here x denotes a ...fire and 
0 a no-fire. On the left, the lines on which there were tests 
are numbered from zero to four. In the table below the chart 
are given the lines numbers, the actual testing currents (in 
amperes), the natural logarithms of the testing; heights, the 
number of explosions, and the number of non-explosions at 
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The following is a chart showing the result of an actual 
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each level. The set of numbers is the same in both columns 
as will always be the case when there are equal numbers of 
zeros and x's . When the total number of zeros is unequal 
to the total number of x's , the number of zeros at a given 
current will not differ by more than one from the number.of 
x's at the next greater current. 

The actual currents were chosen so that their natural 
logarithms were equally spaced, since previous investigations 
of similar explosives revealed that the natural logarithms 
of the critical currents could be assumed to be normally dis- 
tributed. The interval h = .093 in In units was used be- 
cause another experiment with a similar explosive gave this 
value for the standard deviation. 

The mean and standard deviation, will be estimated from 
either the zeros or the x's depending on which occur in 
fewer numbers . In the particular example given in Figure 6, 
either the zeros or the x's may be used, but had there been 
49 zeros and 51 x's, the zeros would have been used. 

Let ng, n 1} n 2 , ••• be the number of zeros (or x's 
as the case may be) on the 0 th , 1 st , 2 nd ••• lines respec- 
tively and let N be the total number of zeros (or x's). 

Let c be the normalized current of the lowest line on which 
there was a test recorded. In the example the x's will 
be used and the zeros disregarded; the numbers just defined 
will then be: N = 50, c - .928, n 0 .= 0, n* = 8, n 2 = 35, 

n 3 - 6, n 4 = .1. 

\ 


18 


The formula for the estimate of the mean Is 


( 1 ) 


m « c + d 


I £ ln i + I 


if the no-flres are used; or 


(<;) 


m 


c + d 




If the fires are used. 


The standard deviation will be determined by computing; 


(3) S = 1.620d< 


N Ji 2 n - (^ir^) 2 
N 2 


+ 0.029 


The sums appearing in (1), (2), and (3) can easily be 
computed on a tabular form: 


i 

n i 

in i 

IS 

0 

0 

0 

0 

1 

8 

8 

8 

2 

35 

70 

140 

3 

6 

18 

54 

4 

1 

4 
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N = 50 

A = 100 

B = 218 


where tne data of the example in Figure 6 have been used, 
and the two sums ^in i and £i 2 n^ have been represented by 

A and B . In terms of A and B the formulas (1), (2j, 
and (3) may be written 
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(5) 

( 6 ) 


in = c + d 


A/N + 


S = l.620d 


NB - A 2 

L N 2 


+ 0.029 


where In (5) the plus sign is to be taken if the no-fires are 
used and the minus sign taken if the fires are used. Equation 
(o) is to be used only if (nb - A 2 )/N 2 > .3 

All computations involving the currents are to be done 
in terms of the normalised currents and only final results 

transformed to actual currents. Thus for the particular ex- 
ample, the mean is 


(7) 


m 

m 


.928 + .073 


100 r 
_ 50 ■ 7 


1.068 


and the standard deviation is 


(8) s 


1.620 (. 093 ) 
.0586 


50(218) - (100) 2 
50 


+ 0.029 


In terms of actual currents, the mean m = 1.068, cor- 
responds to 2.91 amperes; this is not the mean current but 
the median current, that is, the current at which there is 
an even chance that a particular specimen will or will not 
fire. The standard deviation S must always be used in. 
normal units (In units in this case), so there ±r, no point 
in transforming it to amperes. 
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Estimate of the All-Fire and No-Fire Points. 

The estimate of the All-Fire point (A.F.) is given by 
equation (9). 

(9) A.F. - m + 3.09S 

the estimate of the No-Fire point (N.F.) is given by equation 
( 10 ) 

(10) N.F. ■ m - 3*09S 

One-sided confidence limits for A.F. and N.F. 

In order to. calculate the one-sided confidence limits* 
it is necessary to calculate the standard deviation of m 
and the standard deviation of s • These are obtained from 
equation (11) and (12). 

( 13 -) s n ■ (6* ♦ «)/7/r (valid only if d < 3S) 

(12) S 8 ■ (1.1a + .35 a /d)/rtT (valid only if d < 2S) 

the 100a percent one-sided confidence limits for A.F. is 
given by equation (13). 

(13) A.F.C. » A.F. * t #>N-1 ^S* + < 3 . 09 ) 2 Sjj 

where t o ^ N-1 is a value obtained from cumulative t - dis 
tributiori with n-1 degrees of freedom and probability a . 
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The 100a percent one-sided confidence limit for N.P. is given 
by equation (14) 


(14) N.P.C. 


N.P. - t 


i,N-l V S m 


+ ( 3 • 09) 2 




where t is a value obtained from a cumulative t - 

distribution with n-1 degrees of freedom and probability a 


In the example cited above: 


m 


6 (. 0586 ) + .093 
7/50" 


.0090 


= 1-1 (.0586) 4- .3 ( . 0586) 2 /, 093 = 010 

/W 

= 1.296 

1.068 - 3.09 (.0586) = .8864 ; N.P. = 2.42 

amperes 

A.F. = 1.068 + 3.09 (.0586) = 1.2466 ; A.F. = 3.48 


90,49 
N.P. = 


amperes 


Estimate of the No-Pire point with 90 percent confidence 
N.P.C. = 2.31 amperes 

Estimate of the All-Fire point with 90 percent confidence 
A.P.C. * 3.65 amperes 

Maximum Likelihood Solutions 


Since equations (3), (11), and (12) are valid only under 

certain conditions, it is advantageous- to use the maximum 

likelihood method to obtain estimates of m , s , S and 

m 
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S 8 . However, this method requires a detailed computer pro 
gram which is not available at the present time. This pro- 
gram should be available in the near future from: 

The Computation and Analysis Division 
Manned Spacecraft Center 
Houston, Texas 
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